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Chapter  1 
Introduction 


Interfaces  in  polymer  based  carbon  fiber  composites  are  critical  regions  that  are  most  suscepti¬ 
ble  to  delamination  under  static  and  impact  loads.  Mechanical  response  types,  like  compression 
(Fig.  1.1(a)),  bending  (Fig.  1.1(b))  and  impact  (Fig.  1.1(c))  on  layered  composites  are  affected  by 
the  interface  strength  and  toughness,  the  geometry  and  the  extent  of  loading.  Debonding  or  de¬ 
lamination  is  observed  to  be  a  significant  failure  mechanism  in  layered  composites  with  significant 
visible  damage. 


Figure  1.1:  (a)  Compressive  Delamination;  (b)  Bending  Delamination;  (c)  Impact  Delamination 


Design  of  interfaces  is  very  critical  for  layered  materials,  like  fiber  reinforced  composites  and 
bonded  joints,  and  smart  designing  techniques  should  be  developed  to  minimize  the  damage  and 
failure  incurred  by  weak  interfaces.  Improving  the  strength  of  material  is  not  the  only  route  for 
constructing  better  laminated  composites  and  bonded  joints,  instead,  the  interfaces  can  be  engi¬ 
neered  for  specific  types  of  loads  acting  and  the  geometry  of  the  components.  Trial  and  error 
experimental  approach  with  changing  component  geometry  and  load  combinations  will  be  very 
expensive.  Therefore,  a  robust  computational  model  to  drastically  reduce  the  trial  and  error  exper¬ 
iments  for  interface  design  has  to  be  developed.  Interface  strength  and  toughness  can  be  improved 
significantly  by  structural  micro-reinforcements  using  polymer  additive  manufacturing.  A  com¬ 
putational  modeling  approach  is  being  developed  to  obtain  optimized  designs  for  interfaces  and 
upscale  the  influence  of  interface  reinforcements  for  damage  and  failure  resistant  designs,  which 
will  aid  in  improving  the  overall  performance  of  the  composite.  The  strength  of  the  composite 
will  be  increased  with  negligible  increase  in  the  weight  and  cost  associated.  Further,  polymer  ad¬ 
ditive  manufacturing  techniques  will  be  tailored  towards  constructing  these  optimized  microscale 
interface  designs  considering  the  suggestions  obtained  from  the  computational  models. 
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1.1  Research  Objectives 

The  overall  research  objective  of  this  project  are  as  follows: 

1.  Develop  computational  framework  for  designing  optimum  interfaces  for  layered  composites 
to  minimize  damage  and  failure  when  subjected  to  static  and  dynamic  mechanical  loads  as 
well  as  extreme  thermal  loads. 

2.  Develop  a  novel  upscaling  methodology  to  account  for  the  influence  of  interface  stiffeners 
on  the  macroscale  failure  behavior  of  the  composite. 

3.  Investigate  the  feasibility  of  using  additive  manufacturing  for  constructing  micro  scale  in¬ 
terfaces  on  composite  prepregs.  This  will  include  the  construction  of  interfaces  with  spatial 
variation  to  reduce  effects  like,  edge  effects  on  the  interface  delamination  mechanism  and 
damage  due  to  high  intensity  impact. 

4.  Develop  smart  designs  for  bonded  joints  in  aircrafts  using  the  above  developed  computa¬ 
tional  design  methodologies. 

1.2  Background 

Extensive  work  has  been  conducted  by  researchers  that  have  demonstrated  the  influence  of  inter¬ 
faces  on  the  strength  of  the  composite.  Due  to  the  presence  of  fiber  and  matrix  in  a  fiber  reinforced 
composites,  different  failure  mechanisms  characteristic  of  the  constituents  are  observed  when  sub¬ 
jected  to  different  types  of  load  profiles.  For  example,  common  failure  mechanisms  observed 
when  a  laminate  is  subjected  to  compressive  loading  are  kink  band  formation  in  the  lamina  that 
are  aligned  with  the  loading  direction,  matrix  micro-cracking  and  delamination  between  laminae. 
Kink  band  formation,  also  known  as  micro-buckling,  is  an  instability  caused  by  the  axial  loading 
on  the  fibers  in  a  lamina.  In  addition,  a  mismatch  in  the  properties  of  adjacent  layers  causes  stress 
concentration  between  laminae,  resulting  in  the  interface  matrix  rich  region  to  damage  and  create 
new  surfaces  by  releasing  energy.  This  type  of  failure  is  referred  to  as  delamination  failure.  A 
combination  of  both  kinking  and  delamination  is  observed  at  failure.  Lee  and  Waas[l],  Lee  et 
al.  [2],  Vogler  et  al.  [3]and  Pimenta  et  al.[4]and[5],  have  shown  that  micro-buckling  in  the  fibers 
occur  along  with  splitting  when  multi-directional  laminates  are  subjected  to  compressive  loading. 
Lee  and  Waas  [1]  studied  the  effect  of  fiber  volume  fraction  on  the  compressive  failure  mode,  while 
Yerramalli  and  Waas  [6] studied  the  effect  of  fiber  type  and  load  multi- axiality  on  failure.  In  both 
of  these  studies,  energy  released  by  splitting  (debonding)  failure  in  combination  with  kinking  was 
identified  as  contributing  to  the  failure  mechanism.  The  influence  of  interface  delamination  on 
the  compressive  strength  in  multidirectional  fiber  reinforced  laminates  was  demonstrated  through 
validated  computational  modeling  by  the  PI  (Prabhakar  and  Waas  [7]),  where  the  mode-II  shear 
strength  of  the  interface  appears  to  influence  the  compressive  strength  of  the  laminate.  Also,  the 
free  edge  effects  contributed  significantly  to  the  de lamination  response  of  the  composite. 
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1.3  Research  Plan 


The  proposed  work  deals  with  developing  high  fidelity  computational  approach  for  determining 
critical  interfaces  and  smart  design  of  the  same  for  polymer  matrix  composites  and  bonded  joints. 
Also,  a  pioneering  idea  towards  the  building  the  designs  using  additive  manufacturing  is  investi¬ 
gated.  The  steps  mentioned  below  were  proposed: 

•  Task  1:  Develop  interactive  computational  framework  for  designing  structural  interfaces 
for  layered  composites:  Reduced-order  Critical  Interface  Analysis  (RoCIA)  and  Microscale 
Optimized  Design  of  Interface  Truss  Structures  (MODITS) 

•  Task  2:  Develop  a  multiscale  modeling  approach  to  upscale  from  microscale  to  macroscale 

•  Task  3:  Investigate  the  feasibility  of  constructing  structural  interfaces  using  additive  manu¬ 
facturing:  Fused  Deposition  Modeling  (FDM)  type  of  printing  on  Prepregs 

•  Task  4:  Engineering  the  interfaces  using  the  above  formulations  for  aircraft  bonded  joints 

This  is  a  final  technical  report  for  the  project  with  start  date  of  May  1st,  2015  to  August  31st, 
2016.  During  this  period,  the  project  has  mainly  focused  on  Tasks  1  and  4.  Chapters  2  and  3 
describe  the  work  conducted  and  corresponding  findings  within  these  tasks. 

1.3.1  Task  1:  Develop  interactive  computational  framework  for  designing 
structural  interfaces  for  layered  composites 

Under  this  task,  a  computational  approach  for  the  microscale  structural  designs  for  enhanced  inter¬ 
faces  between  layered  fiber  reinforced  composites  is  being  developed.  The  first  step  towards  that 
will  be  formulating  a  reduced  order  finite  element  analysis  (FEA)  approach  to  determine  critical 
interfaces  within  the  composite  when  subjected  to  different  types  of  loads.  Using  the  stress-strain 
information  at  the  interfaces  determined  this  step,  the  optimized  microscale  structural  designs  for 
the  interfaces  will  be  determined  using  a  combination  of  FEA  and  genetic  algorithm  (GA)  ap¬ 
proaches. 

A  novel  modeling  approach,  called  reduced-order  critical  interface  analysis  (RoCIA),  is  being 
developed  to  determine  the  critical  interface  within  layered  composites  to  minimize  the  computa¬ 
tional  cost  associated  with  large  models.  The  framework  will  also  include  the  effects  of  spatial  and 
temporal  variation  of  the  interface  properties  to  determine  the  optimum  interfaces. 

1.3.2  Task  4:  Engineering  the  interfaces  using  the  above  formulations  for 
aircraft  bonded  joints 

The  goal  is  this  task  is  to  first  investigate  the  feasibility  of  incorporating  designs  at  the  bonded 
regions  to  improve  their  bond  strength  using  Fused  Deposition  Modeling  (FDM),  which  is  an 
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additive  manufacturing  technique  for  polymers.  Single  lap  shear  tests  have  been  investigated  to 
explore  the  influence  of  structural  interfaces  at  the  bond  regions  of  single  lap  joints.  The  result 
thus  far  have  been  promising  with  the  apparent  shear  strength  increasing  quadruple  times.  The 
design  incorporated  at  the  bond  regions  on  the  adherents  have  been  analyzed  using  computational 
models  to  study  their  effects  on  the  stress  distribution. 


1.4  Milestones  Achieved:  5/1/2015-8/31/2016 

The  milestones  achieved  so  far  within  this  project  are  as  follows: 

•  Chapter  2:  Reduced-order  Critical  Interface  Analysis  (RoCIA)  to  identify  critical  interlami¬ 
nar  regions  in  multidirectional  laminates  subjected  to  combined  tensile-temperature  loading 
has  been  developed  and  validated. 

•  Chapter  3:  The  feasibility  of  constructing  structural  interfaces  using  additive  manufacturing: 
Fused  Deposition  Modeling  (FDM)  on  laminates  to  improve  their  apparent  shear  strength 
has  been  successfully  demonstrated  through  single  lap  shear  tests  on  bonded  joints 

1.5  Publications 


Journal 


•  M.  S.  Islam,  R  Prabhakar,  “Modeling  Framework  for  Free  Edge  Effects  in  Laminates  under 
Thermo-Mechanical  Loading”.  (Under  Review) 

•  M.  S.  Islam,  P.  Prabhakar,  “Interface  Strengthening  of  Multidirectional  Laminates  using 
Additive  Manufacturing”.  (Under  Review) 

•  R.  Garcia,  E.  Acuna  and  P.  Prabhakar,  “Additive  Manufacturing  for  Improving  the  Strength 
of  Adhesively  Bonded  Single  Lap  Joints  for  Composites”.  (In  Prep.) 

Conference 


•  Md  S.  Islam  and  P.  Prabhakar,  “Interlaminar  Strengthening  of  Multidirectional  Laminates 
Using  Additive  Manufacturing”,  Additive  Manufacturing  of  Composites  and  Complex  -  Ma¬ 
terials,  Materials  Science  &  Technology  Conference  (MS&T16),  Salt  Lake  Cuty,  UT,  Oct. 
23-27,  2016.  (Accepted  for  presentation) 

•  Md  S.  Islam  and  P.  Prabhakar,  “Lree  Edge  Effect  in  Multi-directional  Laminate  Under  Multi- 
axial  Loading”,  American  Society  for  Composites  31st  Technical  Conference  and  ASTM 
Committee  D30  Meeting ,  Williamsburg,  Virginia,  September  19-22,  2016.  (Accepted  for 
presentation) 
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R.  Garcia,  E.  Acuna  and  R  Prabhakar,  “Additive  Manufacturing  For  Bonded  Composite 
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D30  Meeting,  Williamsburg,  Virginia,  September  19-22,  2016.  (Accepted  for  presentation) 
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Chapter  2 

RoCIA:  Novel  Model  for  Free  Edge  Effects 
in  Laminates  under  Thermo-Mechanical 
Loading 

2.1  Introduction 

Interlaminar  regions  in  polymer  based  layered  carbon  fiber  composites  are  the  critical  regions  that 
are  most  susceptible  to  delamination  under  static  and  dynamic  loading.  Mechanical  response  to 
load  types  like  edge-wise  and  through-thickness  compression,  bending,  and  dynamic  impact  on 
layered  composites  is  influenced  by  interlaminar  strength  and  toughness,  the  geometry  and  the  ex¬ 
tent  of  loading.  Debonding  or  delamination  is  observed  to  be  a  significant  failure  mechanism  in 
layered  composites  with  considerable  visible  damage.  Significant  localized  interlaminar  stresses 
occur  at  the  free  edges  of  laminated  composite  materials  due  to  mismatch  in  property  between 
plies,  which  is  referred  to  as  the  “free  edge  effect”  [8].  Accurate  determination  of  stress  distri¬ 
bution  near  the  free  edges  is  very  important  due  to  their  significant  impact  on  delamination  or 
transverse  cracking.  Stress  state  near  the  free  edge  is  three  dimensional  in  nature  and  classical 
lamination  theory  (CLT)  is  unable  to  determine  these  stresses  accurately  [8,  9].  Therefore,  various 
analytical  and  numerical  approaches  such  as  finite  difference,  finite  element,  closed  form  analyt¬ 
ical  approach,  boundary  layer  theories  and  layer-wise  theories  have  been  used  to  calculate  the 
interlaminar  stresses  near  the  free  edges. 

Puppo  and  Evensen  [10]  proposed  the  first  analytical  method  to  determine  the  interlaminar 
stresses  in  a  composite  laminate.  Pagano  [11]  developed  a  higher  order  plate  theory  to  evaluate 
the  interlaminar  stresses.  Few  other  analytical  methods  are:  Perturbation  method  by  Hsu  and  Her- 
akovich  [12];  Boundary  layer  theory  by  Tang  [13],  Davet  and  Destuynder  [14],  and  Lin  and  Ko 
[15];  An  approximate  elasticity  solution  by  Pipes  and  Pagano  [16];  Variational  principle  by  Pagano 
[17],  Lekhnitskii  and  Fern  [18];  Stress  potential  method  by  Wang  and  Choi  [19,  20,  21,  22],  Kas- 
sapoglou  and  Lagace  [23,  24]  and  Webber  and  Morton  [25].  Becker  [26]  used  a  particular  warp 
deformation  mode  to  investigate  the  free  edge  effects  in  laminates.  Kassapoglou  [27]  generalized 
the  principle  of  minimum  complementary  energy  and  force  balanced  method  to  analyze  laminates 
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under  in-plane  and  out-of-plane  loading.  Yin  [28,  29]  used  a  variational  method  using  Lekhnitskiis 
stress  function  for  laminates  under  combinations  of  mechanical  loading.  Later,  Yin  [30]  used  the 
principle  of  complementary  energy  based  on  polynomial  stress  functions  to  evaluate  the  interlam¬ 
inar  stresses  in  laminates  subjected  non-uniform  thermal  loading. 

Pipes  and  Pagano  [31]  developed  the  first  numerical  method  to  solve  the  two  dimensional  gov¬ 
erning  elasticity  equations,  where  finite  difference  method  was  used  to  calculate  the  interlaminar 
stresses  of  long  symmetric  laminate  under  uniform  axial  strain.  Later,  Atlus  et  al.  [32],  Bhaskar  et 
al.  [33]  and  Salamon  [34]  used  three  dimensional  finite  difference  method  to  determine  the  inter¬ 
laminar  stresses  in  angle-ply  laminate.  Wang  and  Crossman  [35,  36]  studied  symmetric  composite 
laminate  subjected  to  uniform  axial  strain  and  thermal  loading  using  finite  element  method.  Her- 
akovich  et  al.  [37],  Isakson  and  Levy  [38],  Rybicky  [39],  Kim  and  Hong  [40],  Icardi  and  Bertetto 
[41],  Lessard  et  al.  [42]  and  Yi  and  Hilton  [43]  also  used  finite  element  method  to  study  free 
edge  effects  in  laminates.  Spilker  and  Chou  [44]  used  hybrid  stress  finite  element  method,  Lee 
and  Chen  [45]  used  layer  reduction  technique,  Robbins  and  Redy  [46]  used  displacement  based 
variable  kinematic  global  local  finite  element  method,  Gaudenzi  et  al.  [47]  used  three  dimensional 
multilayer  higher  order  finite  element  method  to  study  similar  problems. 

From  previous  studies  mentioned  above,  it  is  well  established  that  free  edge  effects  are  dom¬ 
inant  in  multidirectional  laminates  and  lead  to  very  high  interlaminar  stresses  that  prematurely 
initiate  inter-layer  delamination.  Therefore,  determining  interfaces  with  very  large  interlaminar 
stresses  is  critical  for  the  assessment  of  delamination  driven  failure.  By  doing  so,  the  interfaces 
most  susceptible  to  delamination  can  be  determined  and  strengthened  accordingly  during  manu¬ 
facturing  to  reduce  their  susceptibility  to  failure.  Towards  that,  a  Quasi-2D  formulation  within 
the  finite  element  method  (FEM)  framework  is  established  in  this  paper  to  determine  delamination 
prone  interlaminar  regions  in  multidirectional  laminates  subjected  to  thermo-mechanical  loading. 

In  the  current  paper,  the  variational  formulation  presented  by  Martin  et  al.  [48]  is  extended  for 
combined  thermal  and  axial  loading,  referred  to  as  “Quasi-2D”  (Q-2D)  formulation.  The  Q-2D 
formulation  is  implemented  within  the  FEM  framework  to  accurately  determine  the  stress  distribu¬ 
tion  near  the  free  edges  for  cross-ply  ([0/90]s)  and  quasi-isotropic  ([45/  —  45/90/0]s)  laminates. 
Comparison  of  the  current  results  with  previously  published  results  and  corresponding  3D  model 
results  establishes  the  validity  and  accuracy  of  the  Q-2D  model.  This  paper  is  divided  into  the 
following  sections:  Section  2.2  describes  the  mathematical  formulation  for  the  Q-2D  followed  by 
the  details  of  the  implementation  within  the  finite  element  framework  in  Section  2.3.  Discussion 
of  results  from  several  case  studies  are  presented  in  Section  2.4  followed  by  conclusions. 


2.2  Mathematical  Formulation 

Consider  a  laminate  of  length  2 L,  width  2b  and  lamina/layer  thickness  of  h  as  shown  in  Fig.  2.1(a). 
Tensile  load  is  applied  on  the  edges  at  Y1+l  and  along  the  X\  direction,  and  the  edges  at 
So  and  S2 b  are  the  free  edges.  Further,  a  temperature  change  of  AT  occurs  uniformly  within  the 
laminate.  Cross-section  of  the  laminate  consisting  of  N  layers  is  shown  in  Fig.  2.1(b). 

The  stress  components  are  assumed  to  be  independent  of  X\  in  regions  sufficiently  far  from  the 
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dSn-* 


Layer  p+T~ 


Layer  p 


X3  =  +Nh/2 


X,=  0 


x,=  2b 


x3  =  -Nh/2 


(a)  (b) 

Figure  2.1:  (a)3D  Laminate;  (b)  Cross-section  of  a  3D  laminate 


loading  surface  [31,  48],  such  that  the  displacement  field  {U}  can  be  defined  as, 

U1(x1,X2,X3)  =  Ui(x2,X3)  +  6uXi 
U2(x  !,X2,X3)  =  U2(x  2,x3) 

U3(x  i,x2,x3)  =  U3(x2,x3) 


(2.1) 


Here,  en  is  the  uniform  strain  applied  to  the  laminate  along  the  x\  direction.  The  displacement 
field  {U}  and  corresponding  stress  field  {a}  adhere  to  the  following  equations: 


da, 


dxj 


=  0,  Gij  =  aijkhe™h ,  Vz  =  1,2,3  within  H, 


/rr,  l.dUi  dUu  ^ 
f«(C,)=2(^+&-)wlthmn; 

[Ui]  =  0,  [aijTij]  =  0  on  the  interfaces  Yp\  [a^rij]  =  0  on  ^  and  ^ 


2b 


[aijnj\ 


Ft  on  and 

+L  jv- 


%3nj\ 


-Fi  on  and 

-L  - at| 


(2.2) 

(2.3) 

(2.4) 

(2.5) 


Defined  by  equation  (2.1),  introducing  a  trial  field  V  with  en  =  0  and  averaging  equation  (2.2) 
gives, 


da. 


v 


In  dxj 


Vidn  =  0,  VC*  with  i=l,  2,  3. 


Applying  divergence  theorem  to  the  above  equation  yields, 

dVi 
'n  dxj 


aij-^dn=  /  FiVidS~  I  FiVidS 
1  •'H+L 


(2.6) 


(2.7) 


E-i 


Since,  V  does  not  depend  on  x\ : 


BV 

vuxx^dVL  =  0,  VC*  with  l  =  2,  3. 
dxi 


(2.8) 
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Using  al3  =  aijkhe^h  yields, 

P  f)TT  ot  / 

/  ailkh(-^  -  akhAT)—^dfl  =  0,  VU8  for  i,  k,  h=l,  2,  3;  1=2,  3  (2.9) 

Jn  dxh  oxi 

where,  {em}  =  {e}  —  {a} AT. 

Substituting  the  displacement  field  assumption  given  in  equation  (2.1), 

[  aukh(^^  -  akhAT)^Adx1dx2dx3  + eu  [  aim^Adx1dx2dx3  =  0  (2.10) 

Jn  oxh  oxi  Jq  oxi 

Dividing  the  volume  integral  into  two  integrals,  one  along  the  x\  direction,  and  the  other  in  the 
plane  of  the  cross  section  (x2,x3)  gives, 

r+L  r 


I  dx i  [  ~  akpAT)^Adx2dx3- 

J-L  JS  OXp  OX  i 

+L  r  gyi  f+L  f  dVi 

dx i  /  aukiaklAT —1dx2dx3  +  eu  /  dx i  /  aim  — 1dx2dx3  =  0, 

oxi  J_L  Jo  dxi 


(2.11) 


J-L  J  S 

with  1,(3  =  2,  3. 


[  Uiikpi — —aki3AT)-—dx2dx3  —  AT  f  aukiaki— —  dx2dx3  —  en  [  ami— -dx2dx3  (2.12) 
Js  oxp  dxi  Js  dxi  Js  dxi 

Applying  divergence  theorem  on  the  right  hand  side 

[  aukp{^ -  -  akpAT)^Adx2dx3  =  AT  [  aukiotkiViTiids  -  en  [  aimVinids  (2.13) 

Jsp  dx0  dxi  Jds  Jas 

where,  “s”  represents  a  coordinate  that  denotes  the  boundary  dS,  starting  at  the  origin  of  the 
x2  —  x3  axes  for  the  region  S  and  traversing  in  counter  clockwise  direction.  Therefore,  “s”  is  either 
ux2 ”  or  “2:3”  depending  on  the  edge  on  the  boundary  being  traversed. 

The  above  equation  is  modified  to  account  for  layers  with  different  orientation  in  a  multidirec¬ 
tional  laminate  as  given  below, 


A  r  afjp  f)vp  /'  J-L  r 

alk^^-^kp^T)-*Jrdxidx  3  =  ATJ2  aUkiaPkiVfn^ds-enJ2  aunVinIids 

p=1  JSP  °XP  axl  p=1  JdSP  p=1  JdSP 


(2.14) 


Developing  the  right  hand  side  of  2.14,  contribution  of  each  edge 

N  /  pNh/2+ph  f‘Nh/2-\-ph 

ArE  /  a?2fcl  apklVp(2b,x3)dx3-  I 

aVi2kl  a?klV?(0,x3)dx3 

p=i 


I -Nh/2+(j>-l)h 


I -Nh/2+(p-l)h 


N  /  pNh/2+ph  i'Nh/2+ph 

~  Cll  E  4n Vp(2b,x3)dx3-  I  ap211Vp(0,x3)dx3 

p= 1 


(2.15) 


>-Nh/2+(p-l)h 


’  -Nh/2+(p-l)h 
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Contribution  of  each  interface, 


26  , 

Kst1! alV  -  N -  +  ph)dx2 


p= 1  v  70 

7V_1  /  h  \ 

-en  J2y~Jo  (a«ii  _  atfi1)1'f  (^2,  +  Ph)dx2  J 


(2.16) 


and  contribution  of  top  and  bottom 


f  r2b  h  r2b  h  \ 

AT\  Jo  ali3kialklV^X2,~N2^dX2  +  J0  a^iakiVi^(x2:N-)dx2J 

/  r2b  h  r2b  r  h  \ 

-Cr  y~  Jo  ai3ii  Vi1(x2,-N-)dx2  +  a%uV^(x2,N-)dx2J 


(2.17) 


Equation  2.14  is  a  generalized  2D  formulation  which  has  the  displacement  field  along  the 
Xi,  x2  and  x3  directions.  The  input  to  the  above  formulation  is  the  fourth  order  elasticity  tensor  of 
each  layer  of  the  laminate  for  a  linear  elastic  material,  coefficient  of  thermal  expansions,  applied 
external  strain  and  change  in  temperature.  The  effective  loads  calculated  for  a  laminate  are  applied 
to  the  2D  generalized  representation  of  the  laminate  in  the  finite  element  method  model  explained 
in  the  following  section. 


2.3  Implementation  of  the  Quasi-2D  Formulation 

The  above  formulation  can  be  implemented  in  several  ways  within  finite  element  method.  The 
method  developed  in  this  paper  is  to  modify  a  thin  slice  of  a  3D  model  to  behave  like  a  generalized 
2D  model.  A  3D  model  with  a  small  thickness  (1  mm  in  this  case)  in  the  :/'i -direction  is  considered 
as  shown  in  Fig.  2.2.  Multi-point  constraints  are  applied  between  the  front  and  back  face  of  the 
model  such  that  the  displacement  fields  are  independent  of  Xi -direction.  In  Fig.  2.2(a),  m  is  a  nodal 
point  on  the  front  surface  and  n  is  the  corresponding  nodal  point  on  the  back  surface.  These  two 
points  (m  and  n)  are  tied  together  such  that  the  displacements  u,v  and  w  are  identical  for  the  pair 
nodes.  Similarly,  all  the  nodal  points  on  the  front  surface  are  tied  with  their  corresponding  points 
on  the  back  surface,  which  satisfies  the  requirement  of  displacement  fields  being  independent  of  the 
x\ -direction.  This  simulates  the  left  hand  side  of  equation  (2.14).  The  external  loads  established 
by  the  right  hand  side  of  equation  (2.14)  are  applied  to  the  model  on  the  edges  in  the  x2  —  x3  plane 
and  the  interfaces  between  the  layers. 

Eight  noded  linear  hexahedral  elements  (C3D8  in  ABAQUS  notation)  are  used  to  mesh  the 
Q-2D  model  as  shown  in  Fig.  2.3(a).  Very  fine  mesh  is  used  near  the  interfaces  and  free  edge 
to  sufficiently  capture  the  stress  singularities  within  these  regions.  However,  only  one  element  is 
used  along  the  X\  direction  to  enforce  multi-point  constraint  described  above.  Due  to  significantly 
large  number  of  nodes  on  the  front  and  back  faces,  enforcing  multi-point  constraints  manually  is 
cumbersome. 
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(a) 


(b) 


Figure  2.2:  3D  slice  of  a  laminate 


Back  Face 


(a) 


(b) 


Figure  2.3:  (a)  Typical  mesh  in  the  Q-2D  model;  (b)  Working  principle  of  python™  script 


A  python™  script  is  written  that  reads  the  ABAQUS  input  file  and  writes  the  multi-point  con¬ 
straints.  Working  procedure  of  the  python™  script  is  shown  in  Fig.  2.3(b).  The  input  file  obtained 
from  ABAQUS  without  constraints  is  read  by  the  python™  script,  which  identifies  the  correspond¬ 
ing  nodes  from  the  front  and  back  faces.  The  script  then  writes  the  constraint  equations  for  every 
node  in  a  different  file  and  finally  it  rewrites  the  main  input  file,  which  is  capable  of  accessing  the 
equation  file  generated.  Fig.  2.4(a)  shows  a  snapshot  of  the  re-written  input  file  which  includes  the 
equation  file  as  an  “input”  and  Fig.  2.4(b)  shows  a  snapshot  of  the  equation  file  where  the  equation 
constraints  are  written. 

As  an  example,  Fig.  2.5  shows  the  boundary  conditions  for  [45/— 45/90 /0]s  laminate  subjected 
to  an  uniform  axial  strain,  which  is  derived  from  Equation  2.14.  As  mentioned  in  the  previous 
section,  the  input  to  the  Q-2D  model  is  the  4th  order  elasticity  tensor  of  the  individual  layer  and 
the  applied  strain.  The  effective  loads  that  depend  on  the  elastic  constants  are  applied  as  pressure 
along  the  normal  directions  and  surface  traction  in  the  parallel  directions  as  shown  in  Fig.  2.5.  x\, 
X2  and  xt  indicate  the  directions  of  load  being  applied,  and  are  represented  by  the  arrows  showing 
the  directions  of  load  acting  on  different  layers.  The  loads  corresponding  to  xt,  direction  are  not 
shown  in  the  figure  since  the  associated  elastic  constants  are  equal  to  zero  in  the  current  example. 
Accordingly,  the  effective  interface  loads  are  applied  along  the  interfaces  between  the  layers. 
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Figure  2.4:  Snapshot  of  (a)  input  file  and  (b)  equation  file 


eii[(ai3n  -  «mi)4  +  (4m  -  «4n)4  +  (43h 


-45  i— > 
“3311)^3 


Figure  2.5:  Boundary  conditions  for  uniaxial  extension 


2.4  Results  and  Discussions 


The  Q-2D  model  is  applied  to  several  laminates  subjected  to  three  main  types  of  loading:  uniform 
axial  extension  (Section  2.4.1),  uniform  thermal  loading  (Section  2.4.2)  and  combined  axial  and 
thermal  loading  (Section  2.4.3).  The  results  are  compared  with  published  (where  available)  and 
3D  model  results  to  establish  the  reliability  of  the  current  modeling  approach. 


2.4.1  Case  1:  Axial  Loading 

In  this  section,  [0/90]  s  and  [45/  —  45/90/0]  s  laminate  subjected  to  an  uniform  axial  (x-direction) 
strain  is  considered.  The  following  material  properties  were  used:  Eu  =  137.9  GPa,  E2 2  =  E33  = 
14.48  GPa,  G12  =  G 13  =  G23  =  5.86  GPa  and  u12  =  ^13  =  v 23  =  0.21  ,  which  correspond 
to  a  high  modulus  graphite/epoxy  laminate  [49].  The  [0/90] s  cross-ply  laminate  subjected  to  an 
uniform  axial  (,X] -direction)  strain  has  been  previously  investigated  by  Zhang  et  al.  [8],  Tahani  and 
Nosier  [9],  Wang  and  Crossman  [36],  Nguyen  and  Caron  [49],  Carreira  et  al.  [50]  and  Zhen  et  al. 
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[51].  Fig.  2.6  shows  the  comparison  interlaminar  stress  distribution  at  the  0/90  interface  of  [0/90] s 
laminate  subjected  to  a  uniform  axial  strain  with  few  previously  published  results.  Even  though, 
Fig.  2.6(a)  shows  that  the  present  model  predicts  a  slightly  higher  interlaminar  normal  stress  at  the 
edges  compared  to  previous  studies,  the  behavior  is  identical.  The  shear  stress  (023)  distribution  is 
comparable  to  other  models  (Fig.  2.6(b)). 


Figure  2.6:  Distribution  of  interlaminar  stresses  along  the  0/90  interface  of  [0/90] s  laminate  sub¬ 
jected  to  uniform  axial  extension  (a)  normal  stress  <r33  and  (b)  shear  stress  cr23 


For  [45/— 45/90/0]s  laminate,  seldom  results  were  found  in  the  literature  to  compare  the  results 
from  the  current  Q-2D  model.  In  this  case,  a  full  3D  model  was  analyzed  and  the  results  from  the 
Q-2D  model  was  compared  with  the  corresponding  3D  model.  Fig.  2.7  shows  the  comparison  of 
Q-2D  and  3D  model  at  different  interfaces,  which  shows  an  excellent  agreement  between  these 
two  models. 


2.4.2  Case  2:  Thermal  Loading 

In  this  section,  [0/90]  s  and  [45/  —  45/90/0]  3  laminates  subjected  to  an  uniform  temperature  change 
AT  throughout  the  modeling  domain  is  considered.  For  direct  comparison  with  the  previously 
published  results,  the  following  material  properties  were  used:  Eu  =  137.9  GPa,  T22  =  T33  = 
14.48  GPa,  G12  =  G13  =  G23  =  5.86  GPa,  u12  =  u13  =  v23  =  0.21,  =  0.36  x  10"6  °C~l 

and  a2  =  a3  =  28.8  x  10~6  °C~1.  Fig.  2.8  shows  the  distribution  of  cr33  (Fig.  2.8(a))  and  <r23 
(Fig.  2.8(b))  for  AT  =  1  °C.  It  is  observed  from  Fig.  2.8  that  the  current  results  are  in  excellent 
agreement  with  those  determined  by  previous  investigators  [9,  30,  49]. 

Fig.  2.9  shows  the  distribution  of  interlaminar  stresses  for  [45/  —  45/90/0] s  laminate  along 
45/-45  (Fig.  2.9(a)),  -45/90  (Fig.  2.9(b))  and  90/0  interfaces  (Fig.  2.9(c))  subjected  to  an  uniform 
temperature  change  AT  =  1  °C.  It  is  observed  from  Fig.  2.9  that  the  current  results  are  in  good 
agreement  with  the  corresponding  3D  model  results. 
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Figure  2.7:  Distribution  of  interlaminar  stress  of  [45/  —  45 / 90 /0]s  laminate  subjected  to  uniform 
axial  extension  along  (a)  45/-45,  (b)  -45/90  and  (c)  90/0  interface 


Figure  2.8:  Distribution  of  interlaminar  stresses  along  the  0/90  interface  of  [0/90] s  laminate  due  to 
a  temperature  change  AT  =  1  °C  (a)  normal  stress  cr33  and  (b)  shear  stress  023 


2.4.3  Case  3:  Combined  Axial  and  Thermal  Loading 

In  this  section,  [0/90]  s  and  [45/  —  45/90/0]  s  laminates  subjected  to  a  uniform  axial  strain  (en  = 
0.01)  and  uniform  temperature  change  (AT  =  C)  is  considered.  The  material  properties  were 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


(a)  (b) 


(c) 


Figure  2.9:  Distribution  of  interlaminar  stress  of  [45/  —  45/90/0] s  laminate  due  to  a  temperature 
change  AT  =  25  °C  (a)  45/-45,  (b)  -45/90  and  (c)  90/0  interface 

the  same  as  mentioned  in  Section  2.4.2.  Fig.  2.10  and  Fig.  2.11  show  the  distribution  of  interlami¬ 
nar  stresses  for  [0/90] s  and  [45/  —  45/90 /0]  s  laminates  subjected  to  en  =  0.01  and  AT  =  25  °C. 
These  results  are  compared  with  their  corresponding  3D  model  and  is  observed  (Fig.  2.10  and  Fig. 
2.11)  that  the  current  results  are  in  excellent  agreement  with  the  3D  results. 

2.4.4  Efficiency  of  Q-2D  Model 

The  Q-2D  modeling  framework  developed  in  this  paper  utilizes  a  fairly  straight  forward  derivation 
and  an  easy  implement  within  finite  element  method.  It  is  also  highly  efficient  in  computational 
cost  compared  to  the  3D  counterpart  models.  Table  2.1  shows  a  comparison  of  the  computational 
time  required  by  the  Q-2D  and  3D  models  of  [45/  —  45/90/0]  s  laminate  («  same  time  for  all 
loading  cases).  The  total  computational  time  required  by  Q-2D  model  is  approximately  2  minutes 
and  for  the  3D  model  is  approximately  60  minutes  using  the  same  Intel®  Core™  i7-3540M  3.00 
GHz  CPU.  Although,  60  minutes  seems  like  a  reasonable  time  for  a  3D  model,  it  is  for  only  8 
layers  laminate  ([45/  —  45/90/0]  s).  The  time  required  will  scale  up  with  increased  number  of 
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Figure  2.10:  Distribution  of  interlaminar  stress  along  the  0/90  interface  of  [0 / 90]  s  laminate  sub¬ 
jected  to  6n  =  0.01  and  AT  =  25  °C  (a)  normal  stress  cr33  and  (b)  shear  stress  cr23 


Figure  2.11:  Distribution  of  interlaminar  stress  of  [45/  —  45 / 90 /0]s  laminate  subjected  to  en  = 
0.01  and  AT  =  25  °C  along  (a)  45/-45,  (b)  -45/90  and  (c)  90/0  interface 


layers  in  the  case  of  a  thicker  laminate  (for  example  48  layers).  Hence,  the  Q-2D  model  can  save 
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the  computational  time  significantly  while  accurately  predicting  the  interlaminar  stress  distribution. 


Table  2.1:  Comparison  of  Q-2D  and  3D  model 


Model 

No.  of  Elements 

Approximate  Time  (min) 

Q-2D 

23424 

2 

3D 

318786 

60 

2.5  Conclusions 

In  this  study,  free  edge  laminates  subjected  to  uniform  axial  extension  and  uniform  temperature 
change  are  investigated  with  Quasi-2D  model.  This  generalized  plane  strain  model  is  capable  of 
predicting  the  interface  stresses  accurately  compared  to  3D  model.  The  accuracy  and  effectiveness 
of  the  present  Q-2D  model  is  validated  with  the  previously  published  results  and  3D  model.  This 
Q-2D  model  is  found  to  be  capable  of  predicting  the  interlaminar  stresses  for  different  laminates 
subjected  to  uniform  axial  loading  and/or  thermal  loading  in  a  reliable  and  efficient  way. 
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Chapter  3 

Additive  Manufacturing  for  Bonded 
Composite  Joints 

3.1  Introduction 

Automotive  and  aerospace  industries,  amongst  many  others,  are  exploring  routes  to  improve  the 
quality  of  joints  in  an  assembled  component.  Joints  are  critical  regions  in  structures  due  to  the 
stress  concentrations  manifested  compared  to  the  members  of  a  component.  With  increasing  ap¬ 
plications  of  fiber  reinforced  polymer  matrix  composites  (FRPCs)  in  aircrafts,  navy  structures  and 
automobiles,  novel  joining  technology  to  assist  the  fabrication  of  large  components  has  become  a 
priority  for  structural  engineers.  Conventional  materials  such  as  steel  or  aluminum  are  joined  using 
fasteners  and/or  bolted  joints,  which  is  not  favorable  for  FRPCs  as  drilling  or  cutting  of  the  com¬ 
posite  may  damage  fibers  causing  an  adverse  effect  on  their  structural  integrity.  Thus,  adhesively 
bonded  joints  are  becoming  a  viable  option  for  joining  FRPCs. 

Advantages  of  bonded  joints  over  traditional  mechanical  fasteners  are  lower  structural  weight 
and  improved  damage  tolerance.  Despite  these  advantages,  bonded  joints  in  primary  load-bearing 
applications  often  result  in  overdesign  due  to  the  inclusion  of  mechanical  fasteners  for  additional 
safety.  This  is  due  to  the  lack  of  confidence  in  adhesively  bonded  joints  for  composite  joining 
technology.  Mechanics  based  designs  for  bonded  joints  are  necessary  to  facilitate  efficient  use  of 
composites  for  lightweight  applications.  This  paper  presents  the  use  of  additive  manufacturing  to 
improve  the  mechanical  behavior  of  bonded  joints  by  enhancing  the  load  bearing  capacity  of  the 
bond  area. 

Aerospace  industry  was  the  first  to  conduct  work  on  adhesively  bonded  joints  in  the  1970-80s. 
An  extensive  review  on  the  joint  strength  of  adhesive  joints  for  fiber  reinforced  composites  was 
provided  by  Matthews  et  al.  [52].  New  analytical  approaches  and  finite  element  methods  were 
developed  for  analyzing  bonded  joints.  The  bond  strength  for  FRPCs  depends  on  various  parame¬ 
ters,  like,  joint  configuration,  adhesive  properties,  surface  preparation,  test  methods,  environmental 
conditions,  etc.  Of  these,  surface  preparation  is  perhaps  the  most  critical  in  governing  the  qual¬ 
ity  of  the  joint  [53].  According  to  Davis  and  Bond  [54],  surface  treatment  with  favorable  surface 
chemistry  prior  to  bonding  can  result  in  significant  increase  in  the  bond  strength  resulting  in  im- 
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proved  durability  of  adhesive  bonds.  In  order  to  ensure  maximum  bond  strength,  abrasion/solvent 
cleaning  techniques  are  commonly  employed  as  surface  treatment  for  thermoset  composites,  like 
carbon/epoxy  or  glass/epoxy  composites.  Other  viable  options  for  surface  treatments  are,  grit 
blasting,  acid  etching,  laser  treatment,  etc.  A  novel  method  is  proposed  here  that  utilizes  additive 
manufacturing  to  impart  engineered  structure  to  adherent  surfaces  in  order  to  drastically  increase 
the  bond  strength.  It  is  hypothesized  that  mechanics  dictated  gradation  of  interface  modification 
will  further  assist  in  improving  the  strength  and  toughness  of  the  adhesive  bond. 

Bonded  joints  can  have  different  geometries  depending  on  the  type  of  application.  Commonly 
studied  joints  by  previous  researchers  are  single-lap,  double-lap,  scarf,  stepped,  T-shaped,  tubular 
lap,  L-shaped  joints,  etc.  A  detailed  procedure  for  the  design  of  adhesively  bonded  joints,  like 
single,  double  and  step,  for  composites  under  static  and  cyclic  loads  is  given  in  Chamis  and  Murthy 
[55].  Stress  distribution  within  a  joint  dictates  the  strength  of  an  adhesively  bonded  joint,  which 
in  turn  depends  on  the  geometry  of  the  joint  and  material  properties  of  adhesive  and  adherent. 
Therefore,  a  detailed  analysis  of  a  joint  configuration  with  the  corresponding  materials  is  required 
to  evaluate  the  stress  distribution.  The  joint  design  should  be  catered  towards  minimizing  stress 
concentrations  that  catalyze  debonding  type  failure  at  the  composite  joint.  Especially,  peel  and 
cleavage  stresses  should  be  suppressed,  while  maintaining  an  almost  uniform  stress  distribution 
along  the  joint.  The  effects  of  stress  concentrations  are  accentuated  due  to  layered  nature  of  fiber 
reinforced  composites.  The  weak  interfaces  within  the  adherents  can  result  in  their  failure  along 
with  joint  failure.  Therefore,  it  is  essential  to  improve  the  mechanical  properties  of  bonded  joint 
regions  of  a  layered  composite. 

Single-lap  joint  is  one  of  the  most  commonly  used  joint  geometry  due  its  simplicity  and  easi¬ 
ness  to  fabricate.  However,  stress  concentrations  are  manifested  at  the  ends  of  the  bond  overlap  that 
potentially  cause  damage  resulting  in  the  premature  failure  of  a  bonded  joint.  Earlier  researchers 
have  explored  modifying  the  adherent  geometries  (  [56]  [57]  [58]),  like  tapering,  stepping  and 
wavy  lap,  as  possible  options  to  minimize  these  stress  concentrations  to  improve  their  load  bear¬ 
ing  capacity.  However,  a  changing  adherent  shape  is  a  constraint  on  the  component  geometry  and 
may  not  be  favorable  for  fabrication  purposes.  Hence,  there  is  a  critical  need  to  improve  the  bond 
strength  while  maintaining  standard  geometries  of  the  adherent.  Further,  complex  joint  designs, 
like  T-joint  or  Pi-joint,  that  are  used  to  join  composites  require  extensive  investigation  on  the  stress 
distribution  and  its  influence  on  the  overall  structural  integrity  of  the  component. 


Motivation 

Interface  design  is  very  critical  for  bonded  joints  in  FRPCs  and  a  smart  designing  technique  should 
be  developed  to  minimize  the  damage  and  failure  incurred  by  weak  bond  interfaces.  It  is  hypothe¬ 
sized  that  improvement  in  strength  and  toughness  can  be  achieved  by  imparting  structural  texture 
at  the  adhesively  bonded  joint  interfaces.  Hence,  novel  technology  using  additive  manufacturing 
at  the  bond  interfaces  has  been  explored  in  this  paper. 

The  research  reported  in  this  paper  was  split  into  the  following  main  steps:  1)  Conduct  single 
lap  shear  (SLS)  tests  on  pure  epoxy  adhesively  bonded  joints  to  establish  a  baseline  for  appar¬ 
ent  shear  strength;  2)  Analyze  the  bonded  joint  systems  computationally  to  determine  the  stress 
distributions;  3)  Modify  the  bond  design  computationally  to  understand  their  influence  on  the  dis¬ 
tribution  of  stresses;  4)  Enable  the  optimum  designs  at  the  bond  regions  accurately  using  fused 
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deposition  modeling  within  polymer  additive  manufacturing;  5)  Conduct  single  lap  shear  (SLS) 
tests  on  the  new  designs  to  check  the  influence  of  the  new  bond  interface  designs  on  the  apparent 
shear  strength. 

The  outline  of  the  paper  is  as  follows:  The  manufacturing  process  for  pure  adhesively  bonded 
joints  and  additively  manufactured  joints  is  explained  first,  followed  by  the  computational  model¬ 
ing  of  altered  bond  designs.  The  experimental  approach  for  testing  the  single  lap  joints  is  described 
next  followed  by  results,  discussion  and  conclusions. 


3.2  Manufacturing 

Single  lap  shear  adhesively  bonded  joints  were  examined  in  this  paper  following  the  ASTM  D5868- 
01  standard.  The  manufacturing  process  of  the  bonded  joints  can  be  divided  into  two  parts:  (1) 
development  of  bonded  joints  using  pure  adhesive  only,  and  (2)  development  of  the  bonded  joints 
incorporating  3D  printed  designs  on  the  adherents.  The  joint  preparation  process  for  both  cases 
will  be  discussed  next. 


3.2.1  Pure  Adhesive  Joints 

Two  carbon  fiber  laminate  adherents  were  bonded  together  with  LOCTITE  Hysol  E-120HP  epoxy 
that  was  cured  at  room  temperature  for  24  hours.  The  adherent  substrates  used  were  water-jet  cut 
from  a  DragonPlate  solid  carbon  fiber  laminate  (3/31”  Matte  12”x  24”).  A  total  of  5  joints  with 
dimension  shown  in  Fig.  3.1  were  fabricated  and  tested. 


0.76  mm  (0.03  in)  Bondline 


Figure  3.1:  Dimensions  for  single  lap  joint  as  per  the  ASTM  standard. 


Bond  line  thickness  of  0.76  mm  was  achieved  using  glass  microspheres,  which  provided  a 
bond  line  control  of  700-800  microns.  A  mass  ratio  of  10:1  between  adhesive  and  microspheres 
was  implemented.  The  microspheres  were  mixed  thoroughly  with  the  resin  prior  applying  on  to 
the  substrates  over  an  overlap  area  of  6.45  cm2.  Pressure  was  exerted  on  the  overlap  area  using 
c-clamps  (see  Fig.  3.2).  In  order  maintain  the  same  applied  pressure  on  each  joint  fabricated,  it 
was  ensured  to  maintain  the  equal  number  of  threads  above  the  c-clamp  overhead  for  each,  thus, 
regulating  equal  pressures  for  all  samples. 
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Figure  3.2:  C-clamps  were  used  for  applying  constant  pressure. 

3.2.2  3D-Printed  Adhesive  Joints 

Manufacturing  3D-printed  adhesive  (3D-PA)  joints  using  fused  deposition  modeling  (FDM)  re¬ 
quired  printing  the  designs  on  to  the  carbon  fiber  adherents  (or  substrates)  prior  to  fabricating  the 
joints  using  epoxy.  FORTUS  900mc  Stratasys  machine  was  the  printer  used  located  within  the 
W.  M.  Keck  Center  for  3D  Innovation  at  the  University  of  Texas  at  El  Paso.  The  material  used 
to  create  interface  designs  was  ABS-M30,  which  will  be  referred  to  as  “model”  material  in  this 
paper.  ABS-M30  is  up  to  25-70  percent  stronger  than  standard  Stratasys  ABS  material,  possessing 
greater  tensile,  impact,  and  flexural  strength.  The  model  material  was  extruded  at  a  temperature  of 
315  C,  the  support  extrusion  temperature  was  320  C,  and  the  chamber  temperature  was  95  C.  Note 
that  these  temperatures  vary  for  different  model  materials,  such  as  polycarbonate  (PC)  or  ULTEM 
9085.  Prior  to  starting  printing  the  model  material  onto  the  substrates  (or  adherents),  they  were 
exposed  to  the  chamber  temperature  for  ABS-M30  for  fifteen  minutes  to  check  for  any  damage 
imparted  due  to  the  chamber  temperature.  The  substrates  showed  no  degradation  in  the  material 
property  and  the  printing  process  was  initiated  next. 

Due  to  the  lack  of  a  standard  manufacturing  procedure  for  3D  printing  over  carbon-fiber  sub¬ 
strates,  a  new  methodology  was  developed  and  followed  for  printing  above  the  nominal  thickness 
of  the  substrates.  3D  printing  machines  are  generally  calibrated  in  the  x,  y,  and  z  directions  for 
printing  over  their  corresponding  platforms,  but,  lack  the  same  for  composite  substrates. 

An  overview  of  the  fabrication  procedure  is  described  here.  First,  a  design  resembling  teeth- 
shaped  lines  was  developed  using  CAD  software.  The  slice  height  parameter  was  set  to  0.007  in, 
which  derived  in  a  four  layer  stack  up  for  the  desired  0.03  in  bond  control  thickness.  Because  of 
this,  no  microspheres  were  mixed  with  the  resin  since  the  bond  line  control  was  achieved  by  the 
printed  design  itself.  These  lines  are  intended  for  reinforcing  the  through-thickness  interface  of 
the  joint,  resisting  shear  from  any  mode  II  stimuli  applied  -  this  will  be  discussed  later.  After  this, 
an  enclosure  box  with  the  nominal  thickness  (z  direction)  of  the  substrates  was  added  below  the 
teeth-shaped  lines  to  ensure  two  things:  (1)  that  lines  printed  exactly  over  the  enclosed  box  where 
the  substrate  will  be  placed  and  (2)  that  it  constraints  the  movement  of  the  composite  in  the  three 
directions.  This  box,  at  first,  was  filled  with  support  material,  but  was  removed  later  using  Stratasys 
Insight  software  (see  Fig.  3.3).  A  pause  was  inserted  in  the  beginning  layer  of  the  teeth  lines  -  this 
way,  once  the  box  was  built  with  no  support,  the  printer  would  automatically  stop,  allowing  the 
user  to  place  the  substrate  inside  the  box,  and  resume  the  printing  job  at  hand.  Printing  time  of  the 


24 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


final  .stl  file  used  was  approximately  ten  minutes 


Figure  3.3:  (a)  Isometric  view  of  enclosure  box  with  teeth  lines;  (b)  Side  view  zoom  of  area  of 
interest. 


Figure  3.4:  (a)  Isometrc  view  of  3D  printed  output;  (b)  Visible  gaps  in  side  view  of  lap  joint. 

Fig.  3.4  shows  the  side  view  of  the  printed  specimens.  The  spacing  between  the  design  lines 
were  adjusted  such  that  there  is  enough  space  for  a  resin  layer  between  them  to  assist  perfect 
bonding  in  the  lap  joint.  Overall,  the  3D  printing  process  involved  the  following  steps: 

•  Prepare  CAD  file  of  the  design  and  enclosure  box  using  SolidWorks. 

•  Pre-processing  of  CAD:  deleting  support  and  adding  a  pause  layer  using  Insight  software. 

•  Build  the  output  file  in  the  FORTUS  900mc  machine,  and  insert  the  substrate  when  the  pause 
layer  begins.  Then,  resume  build. 

•  If  bond  strength  between  ABS-M30  and  carbon  fiber  is  as  desired,  proceed  to  build  next 
file.  If  not,  calibrate  machine  in  the  z  direction  and  re-build  until  desired  bond  strength  is 
reached. 

Upon  printing  the  designs  on  the  substrates,  the  procedure  for  fabricating  the  single  lap  joints 
was  identical  to  the  pure  adhesive  joints  case,  except  microspheres  were  not  used  here  as  the 
designs  provided  the  required  bond  size. 


3.3  Computational  Approach 

A  computational  model  of  the  bonded  joint  specimen  was  developed  (see  Fig.  3.5)  to  determine 
the  stress  distribution  in  the  bond  regions.  Commercially  available  finite  element  method  software 
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Figure  3.5:  (a)  Undeformed  finite  element  model;  (b)  Deformed  finite  element  model. 


“ABAQUS”  was  used  to  model  the  system.  Due  to  the  complexity  of  the  structure  within  carbon 
fiber  laminates,  a  simplified  homogenized  model  for  the  adherents  was  developed. 

The  effective  mechanical  properties  of  the  carbon  fiber  laminate  (adherents),  adhesive  and  the 
model  material  (3D  printed)  are  shown  in  Table  3.1. 


Material 

Young’s  Modulus  (PSI) 

Poisson’s  Ratio 

Adherents 

30,000,000 

0.3 

Adhesive 

1,500,000 

0.3 

Reinforcement  ABS-30M 

400,000 

0.3 

Table  3.1:  List  of  materials  and  their  assumed  mechanical  properties 

The  first  model  generated  was  the  bonded  joint  specimen  with  pure  epoxy.  The  loading  condi¬ 
tion  applied  was  a  0.04”  displacement  in  tension  at  both  ends  of  the  specimen  along  the  x-direction. 
Maximum  principal  stress  and  strain  were  plotted  along  the  centerline  of  the  bonded  region  as 
shown  in  Fig.  3.6. 


Figure  3.6:  (a)  Maximum  principal  strain  along  center  path;  (b)  Maximum  principal  stress  along 
bonded  region  center  path 

This  was  followed  by  modeling  the  bond  region  with  different  designs  as  shown  in  Fig.  3.7. 
Four  distinct  models  with  different  spacing  between  the  3D  printed  lines  were  developed  and  simu- 
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lated.  In  the  pure  epoxy  model,  the  bonded  regions  experienced  significant  stress  and  strain  towards 
the  edges  of  the  bonded  region  [59].  Due  to  these  significant  stresses  and  strains,  the  material  is 
more  likely  to  initiate  failure  in  these  regions.  The  prime  objective  of  the  designs  imparted  to  the 
bonded  interface  is  to  evenly  redistribute  the  stresses  and  strains  to  effectively  reduce  the  risk  of 
failure,  thus  improving  the  strength  of  the  adhesively  bonded  joints. 


Figure  3.7:  (a)  Model  1;  (b)  Model  2;  (c)  Model  3;  (d)  Model  4. 

The  stress  distribution  was  determined  for  each  case,  and  the  results  are  shown  in  Fig.  3.8. 


Figure  3.8:  (a)  Principal  strain  for  pure  epoxy  and  4  models;  (b)  Principal  stress  for  pure  epoxy 
and  4  models. 


It  is  hypothesized  that  the  models  with  uniform  stress  distribution  and  lower  edge  stresses 
compared  to  the  pure  adhesive  case  are  expected  to  display  higher  resistance  to  bond  failure.  This 
will  be  validated  through  SLS  tests  of  the  printed  samples,  where  it  is  expected  that  higher  shear 
strengths  will  be  observed  for  those  models. 
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3.4  Experimental  Approach 


Single  lap  shear  tests  were  conducted  to  determine  the  apparent  shear  strength  of  the  single  lap 
joint  bonded  specimens.  Pure  epoxy  adhesive  and  3D-PA  joints  were  tested  in  the  same  manner. 

3.4.1  Mode-II  Tensile  Testing  of  the  Lap  Shear  Joints 

Tensile  tests  were  performed  on  an  INSTRON  5696  machine,  assuring  safety  conditions  and  fol¬ 
lowing  protocol  regulations  [60]  and  specifications  mentioned  in  ASTM  D5868-014  [61].  The 
initial  grip  separation  was  75  mm  with  25.4  mm  minimum  grip  length  for  samples  at  each  end. 
The  specimen  loading  rate  was  13  mm/min.  Five  samples  for  each  bond  type  were  prepared  as 
mentioned  before  and  tested.  Fig.  3.9  shows  the  experimental  set  up  for  the  lap  shear  test. 


Figure  3.9:  Experimental  set  up  of  shear  lap  tests 


3.5  Results  and  Discussion 

The  SLS  data  obtained  from  the  INSTRON  machine  was  post-processed  to  determine  the  apparent 
shear  strength  of  the  bonded  joints.  According  to  the  ASTM  standard,  individual  peak  load  values 
(KPa),  failure  type,  and  averages  by  maximum  and  minimum  values  were  determined.  Failure  type 
of  the  lap  joints  was  decided  based  on  the  criteria  shown  in  Fig.  3.10. 

The  results  for  the  pure  epoxy  bonded  joints  are  discussed  first  followed  by  the  3D-PA  bonded 
joints.  A  discussion  of  the  comparison  of  the  outputs  from  all  cases  will  be  conducted  in  a  later 
section. 


3.5.1  Pure  Adhesive  Bonded  Joints 

The  peak  loads  determined  from  the  SLS  tests  for  PA  bonded  joints  are  given  in  Table  3.2.  The 
shear  stress  in  each  case  was  calculated  by  dividing  the  peak  load  by  the  bond  area. 
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Fiber-Tear  Failure  Stock-Break  Failure 

Figure  3.10:  Failure  types  in  lap  shear  tests. 


Joint 

Peak  Load  (N) 

Shear  Stress  (MPa) 

Failure  Type 

1 

680.38 

35.15 

Adhesive 

2 

1256.08 

64.90 

Adhesive 

3 

808.92 

41.79 

Adhesive 

4 

13546 

699.88 

Cohesive 

5 

1130.91 

58.43 

Adhesive 

Table  3.2:  Post-processing  results  for  pure  epoxy  joints. 


Adhesive  and  cohesive  failure  were  the  two  key  failure  modes  observed  from  the  five  joints 
tested.  Adhesive  failures  can  be  observed  in  Fig.  3. 1 1(a),  where  the  bond  region  separates  from  the 
adherent  completely.  Note  that  the  glass  microspheres  are  visible  in  this  figure.  Cohesive  failure 
was  observed  in  one  case  only  as  shown  in  Fig.  3.11(b),  where  the  failure  region  passed  through 
the  bond  material.  Cohesive  failure  was  ignored  due  to  its  rare  occurrence  within  the  5  pure  epoxy 
bonded  joints  tested. 


Figure  3.11:  (a)  Top  view  of  adhesive  failure;  (b)  Top  view  of  cohesive  failure. 


3.5.2  3D-Printed  Adhesive  Bonded  Joints 

The  SLS  test  results  for  the  four  3D-PA  bonded  joints  will  be  discussed  individually.  The  data 
format  is  the  same  as  in  the  case  of  PA  adhesive  joints,  as  per  the  standard. 

3D-PA  Joints:  Model  1 

Table  3.3  represents  the  peak  load,  shear  strength  and  failure  types  for  3D-PA  Model- 1  adhesive 
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joints. 


Joint 

Peak  Load  (N) 

Shear  Stress  (MPa) 

Failure  Type 

1 

4102.82 

212.53 

Adhesive 

2 

5254.18 

272.18 

Cohesive 

3 

3755.47 

194.54 

Adhesive 

4 

3943.09 

204.26 

Adhesive 

5 

4610.36 

238.83 

Adhesive 

Table  3.3:  Post-processing  results  for  model  1  joints. 

Model- 1  3D-PA  bonded  joints  displayed  adhesive  and  cohesive  failures  at  the  bond  region  as 
shown  in  Fig.  3.12.  Note  that  the  3D  designs  are  visible  in  both  types  of  failure.  It  should  be  noted 
that  the  corresponding  peak  loads  and  stresses  increased  significantly  compared  to  the  PA  joints. 


Figure  3.12:  (a)  Top  view  of  adhesive  failure;  (b)  Top  view  of  cohesive  failure. 
3D-PA  Joints:  Model  2 


Joint 

Peak  Load  (N) 

Shear  Stress  (MPa) 

Failure  Type 

1 

6551.08 

338.47 

Cohesive 

2 

4778.79 

246.9 

Adhesive 

3 

6851.07 

353.97 

Cohesive 

4 

7755.10 

400.68 

Cohesive 

Table  3.4:  Post-processing  results  for  model  2  joints. 


3.5.3  Conclusion:  Pure  Adhesive  vs.  3D-Printed  Adhesive  Bonded  Joints 

This  initial  study  has  shown  that  the  structure  imparted  to  the  bond  regions  using  3D  printing 
has  significantly  improved  the  shear  strength  of  the  bonded  joints.  Current  ongoing  work  focuses 
on  printing  the  remaining  designs  that  were  computationally  analyzed  to  test  the  hypothesis  stated 
before.  Further  study  will  investigate  the  influence  of  optimized  design  on  the  strength  of  pi-joints. 
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